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Theoretical models which present flavor changing neutrino interactions and simultaneously prevent 
this particle from acquiring any mass exist. We discuss some of them and their predictions for 
neutrino oscillations in matter which can account for the solar neutrino anomaly as well as the 
zenith-angle dependence of the deficit of atmospheric neutrinos observed by the SuperKamiokande 
experiment without invoking, therefore, neutrino masses nor mixing. 



I. INTRODUCTION 

Flavor changing neutrino interactions (FCNI) can in- 
duce neutrino oscillations in matter. This phenomenon 
was first investigated by Wolfenstein |l| who pointed 
out that interactions in a medium modify the disper- 
sion relations of particles traveling through. Wolfen- 
stein effect generates quantum phases in time evolution 
of phenomenological neutrinos eigenstates which conse- 
quently can oscillate. These oscillations can be res- 
onantly enhanced even if neutrinos are massless and 
no mixing in the vacuum exists [^J^. Similar phe- 
nomenon happens when neutrino masses are introduced 
and mixing angles in matter are induced by flavor chang- 
ing interactions Wolfenstein effect has been in- 
voked to obtain a good fit of the solar neutrino observa- 
tions 1^ when FCNI is assumed in the massless neutrino 
context [pUq] or when the neutrinos are assumed mas- 
sive [^I^J^. Furthermore it was recently pointed out that 
the SuperKamiokande results Q showing a strong zenith- 
angle dependence of the /i flux induced by atmospheric 
/i- neutrinos, which have been interpreted as an indication 
of neutrino oscillations ||9|,p^ , can also be understood as- 
suming non-resonant FCNI with ordinary matter in the 
Earth |ll|,|l|l. 

The presence of flavor changing neutrino-matter inter- 
actions implies a non-trivial structure for the neutrino 
evolution Hamiltonian in matter even if massless neu- 
trinos and no mixing in the vacuum is assumed. The 



evolution equations describing the Va i^/s transitions 
{a, /? = e, /i, r are flavor indices) are given by 0,0|: 



dr 



a/3, 



(1) 



where ^'q^ = vp) ; Va = I'ai'i') and = i^^(r), are 
the probability amplitudes to find these neutrinos at a 
distance r from their creation position, \/2 GFnf{r)e^p 
is the flavor-changing Va + f i^p + f forward scatter- 
ing amplitude with the interacting fermion / (charged 
lepton, c?-like quark or u-like quark) and \/2 GFnf{r)e'^p 
is the difference between the flavor diagonal Va — f and 
Vfj — f elastic forward scattering amplitudes, with nf{r) 
being the number density of the fermions which induce 
such processes. In all cases of practical interest, electron- 
neutrinos will coherently scatter off the electrons present 
in matter through standard electroweak charged currents 
which introduce non trivial contributions to the neutrino 
evolution equations. These contributions are taken into 
account by the term \/2GF5eane(r) in Eq. (|l|). There- 
fore, when the electron neutrino evolves in matter (a — e) 
and the fermion which the neutrino scatter off is not 
electron (/ 7^ e), a resonance can occur if the condition 
('^pTifir) = nf,{r) is fulfilled at least in one layer of the 
matter crossed by the neutrinos ||] . In the case where no 
electron-neutrino participates in the oscillation, no reso- 
nance can happen. 
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In any case, flavor oscillations induced by FCNI can be 
relevant for understanding a considerable range of exper- 
imental data involving neutrinos which cannot be fitted 
by using only standard electroweak inputs. The possi- 
bility that these oscillations occur even in the absence of 
neutrino masses and mixing open a new perspective for 
interpretation of current neutrino data. Our purpose in 
this paper is to investigate how feasible this picture can 
be in several models which allow FCNI. 

II. FCNI: A GENERAL REVIEW 

Although FCNI do not appear in the standard model 
of the electroweak interactions they are quite naturally 
present in many of its extensions. Most of them, how- 
ever, do not naturally admit strictly massless neutrinos. 
An example of these models are the i?-parity broken su- 
persymmetric models [^3| which were invoked in the rel- 
evant analyses for the solar and atmospheric neutrinos 
previously cited P,^,pl| . In these models, the same FCNI 
interactions that induce the e^^ and e'^^ contributions to 
the Hamiltonian in Eq. (|l|) induce also contributions to 
the neutrino mass at the 1-loop level. A fine tuning in 
the model parameters is needed in order to keep neutrino 
masses negligible while endowing e;^^ and e'^p with ap- 
propriate values to fit the solar and atmospheric neutrino 
anomalies. 

There are extensions of the standard SU{2) (g) U{1) 
model or of its SU{5) grand unification extension in 
which neutrinos get naturally a small Dirac mass [ p^ . 
In order to obtain such effect extra symmetry and fields 
are added to the minimal models. For instance, in the 
SU{2) (g) U{1) model besides the right-handed neutrino 
vr, a pair of extra singlets sl and sr and one scalar sin- 
glet are added. The extra symmetry imply that for each 
family there are one light and one heavy neutral leptons. 
In this case the FCNI needed to induce non-zero e-^ „ and 

ctp 

^^'^ mediated by the gauge vector boson W and for 
this reason the only free parameters are the mixing an- 
gles. This has to be confronted with constraints imposed 
by the experimental limits on lepton number violating 
decays such as /x — > 67 in a similar way we will do in 
Sec. 

Also in this context it is worth to mention other models 
like the left-right symmetric one and their SO(IO) grand 
unified extensions . In this case three of the neutri- 
nos can be massless at arbitrary order in perturbation 
theory. However, the size of flavor changing neutral cur- 
rents in the neutrino sector depends on the value of the 
mass scale related to the pattern of symmetry breaking 
that reduces S'O(IO) to SU(Z)c ® SU{2)l (g U{1). Only 
if 50(10) breaks first to SU{2)l «) SU{2)r (g) U{1) it is 
possible to have intermediated mass scale related to the 
energy at which SU{2)r is broken. In general, grand- 



unified theories present both massless neutrinos and fla- 
vor changing interactions . Nevertheless the effective 
flavor changing neutral currents are negligibly small since 
they are inversely proportional to the squared mass of 
the exotic vector boson which is fixed by the condition 
of unification. This is the case of all non-supersymmetric 
grand-unified models since in these models there exist the 
hierarchy problem. 

A version of the Zee's model [0 in which only one ex- 
tra singly charged scalar singlet (besides the usual dou- 
blet) and right-handed neutrinos are added was consid- 
ered in Ref . jl^ . In the simplest case the constraint com- 
ing from the muon decay implies that the atmospheric 
neutrino flux is reduced at most 20%. A way to overcome 
this difficulty is to add another singly charged scalar sin- 
glet. In this case the neutrino masses are also arbitrary 
and the constraints of the muon and other leptonic decays 
can be evaded. Notice that if we do not add right-handed 
neutrinos FCNI take place only for anti-neutrinos. 

Nevertheless there are models where the mechanism 
pointed out in Refs. |^ and can be exactly realized, 
keeping neutrinos strictly massless, while FCNI can exist 
with the required intensity. Basically, vanishing neutrino 
masses are guaranteed in these models due to the con- 
servation of the total leptonic number L. Then, by im- 
posing that there are no right-handed neutrinos to avoid 
Dirac masses, no vacuum expectation values associated 
with neutral scalars to prevent spontaneous violation of 
L and, consequently, no majoron-like Goldstone boson is 
generated, neutrino masses can be kept vanishing. Obvi- 
ously, no L violation in the scalar sector is also allowed, 
like as in the i?-parity violating supersymmetric mod- 
els, since in this case neutrino masses can be generated 
radiatively |l9| ]. 

In this paper we show that these conditions are nat- 
urally fulfilled in some theoretical scenarios and there- 
fore FCNI can coexist with strictly massless neutri- 
nos. In Sec. |lll| we consider the gauge model based on 
SU{?,)c (g 5t/(3)i, g) U{1)n symmetry ( 331 model) |^] 
and in Sec. |l^ we consider the multi-Higgs doublet ex- 
tension of the standard model [ pT[ |. 

III. FCNI IN 331 MODELS 

Although the standard model of electroweak interac- 
tions accommodates all the present experimental results, 
it is not able to give answer for some questions in parti- 
cle physics. One of these questions is the family replica- 
tion problem, which can have an elegant solution in the 
simplest chiral extension of standard model, the gauge 
model based on S'[/(3)c (g SU{3)l (g U{1)n symmetry 
which broken to SU{3)c (g SU{2)l <g U{1)n in some en- 
ergy scale higher than the Fermi one. In the 331 model of 
Ref. neutrinos are kept massless as long as no extra 
right-handed neutrinos are added to the particle content 



of the model and the total Icpton number L is conserved 
in the full Lagrangian [Q. We will be assuming here 
that this is the case. The 331 phenomenology was stud- 
ied in many different contexts but for the first time 
it is applied for neutrino phenomenology. The Yukawa 
interactions in the lepton sector are [ pO|p^ 
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1 
1 

^2 



2 Vl JC'LRlRvt - 2 II IC'll'''rV2 + H.C., (2) 



where /Cll 



^L^^L^ ^RR 



ErGEr, IClr 



IC'^^ = eIg'Er, K,'^^ = eIg'EI- G and G" are sym- 
metric and antisymmetric (they can be complex) Yukawa 
matrices. The symmetric ( antisymmetric) property of 
the matrices above are because the correspondent Higgs 
field is in triplet (sextet) representation 22|. Erx 
are the right- and left-handed mixing unitary matrices in 
the lepton sector relating symmetry eigenstates (primed 
fields) with mass eigenstates (unprimed fields) p4[] : 



t H 
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FIG. 1. Scalar contributions to FCNI in the 331 model. 
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FIG. 2. Scalar contributions to anti-neutrinos FCNI in 331 
model. 



1'l = ElIl, I'r^ErIr, v']^=Elvl, 



(3) 



where we have redefined the neutrino fields so that there 
is no mixing in the current coupled to the gauge vector 
boson W . Note, therefore, that while /Cll and JCrr are 
symmetric matrices and K'j^j^ is an antisymmetric matrix, 
no symmetry relation appears in IClr and /C^^. None of 
the couplings in Eq. (^) depends directly on the charged 
lepton masses and all matrices in Eq. (|^) are not unitary. 
When G and G' are real matrices the matrices JCll and 
ICrr (and the respective primed matrices) are hermitian. 
All scalars in Eq. (|^) are still symmetry eigenstates. 

The interactions of the leptons with the gauge vector 
boson V [E0[ induce also FCNI and are given by 



Cv = iR.riCv'i.V- + H.C., 



(4) 



with /C — [g^/ \/2)e\^E1 being a unitary matrix and 
is the coupling of the gauge vector boson V with the 
leptons. This is not the g coupling of standard model 
gauge vector boson W. 

We can now write the expressions for the parameter 
entering in Eq. (|^) and for the analogous parameter 

—f 

"^hen anti- neutrinos are evolving instead of neutri- 
nos, in the light of this model. We present, in Figs. 1, 2 
and 3 the diagrams which illustrate scalar contributions 
for neutrinos, scalar contributions for anti-neutrinos and 
vector contributions for anti-neutrinos, respectively, in 
the 331 model. From Figs. |l|, | and | and Eqs. @ and 
(^, we write the nonstandard contributions to e^^ and 
^ap (^^^ superscript e denotes that the relevant interac- 
tions are those with electrons): 
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AV2GfM^ 



{^'lR)^i^'LR)piy 



(5a) 



where x = M^/Mfj^ , y = M^/M^ , and 



aV2GfM^ 



w 
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(5b) 



where z' = M^/M^, x' = M^^/M^^, y' = M^/ml^. 
Note also that we are assuming that all the matrices en- 
tering above are real. In case they are complex, only their 
real part will contribute to e^^ and e^^. In this paper 
we always assume the convention that Greek subscript 
stands for neutrinos and Arabic subscript for charged 
leptons. 
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FIG. 3. Vector contributions to anti-neutrinos FCNI in 331 
model. 

From the above Eq. (^^ we can explicitly calculate 



e^r and e^^ 



e'jj^^ that enter Eq. (|1|) for the case 
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of transitions - 
spheric neutrinos: 



which can be relevant to atmo- 



LRI, (^Lfl)rl y 



(6a) 



and 



- (4.) 



1^ 
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(6b) 
transitions 



Similar for anti-neutrinos, where i 
are at play, we obtain from Eq. (^): 

(/C[^)^^ + (/Ci^Ji^ (/C2j^^ y'} , (7a) 

and 

1 



{[(/C)i,(/Ct),i-(/C)i^(/Ct)^i] z' 



2/ 



(7b) 



Notice that neutrinos and anti- neutrinos will not, in 
general, interact in the same way through the boson 
fields. This is not a CP violation effect but rather only 
reflects the fact that these particles have different inter- 
actions with matter in this theoretical scheme. In fact, if 
all couplings and the vacuum expectation values are real, 
CP is conserved in all the vertices of the model. The pure 
gauge boson interactions also conserve CP if the SU{3)l 
gauge and U{1)n vector bosons W"" with a = 1, ...8, and 
B, respectively transform as 

{W'^,Wl Wl, Wl Wl, Wl Wl, B,0 ^ (8) 
- (W'", -W^", -H/*", W^^", -M^"", H/^", W*", B"). (9) 

It means that the physical fields transform as 

CP 



{w+,v+,u;t+,A,.,z^,z'^) 

(W''',-V''~,-C/''-~,^f,Zf ,Z'f). 



(10) 

(11) 



where W is the standard model gauge vector boson, V 
is a gauge vector boson , J7 is a doubly charged gauge 
boson and and Z'^ are neutral gauge vector bosons. 

In order to obtain the relevant FCNI parameters for 
solar neutrino analysis, i.e., transitions involving and 



i^r, it is sufficient to change appropriately the indexes 
/I — > e in the Eqs. ( |6al ) and (6b). Notice also that in 

I —/- I F-^ I o n A OTt-nil o T-lTr I ^'-^ I — /- I Z^-^ I \zi 



general |e^^| ^ |e^g|, and similarly \e^g\ ^ \e 



e^p I , we need both 

the contributions of the gauge vector boson V to be small 
and a fine tuning among the Yukawa couplings and the 
mixing angles. 

Another important remark concerning the model is 
the following. In order to get to all charged leptons 
an appropriate mass only the sextet of scalars is nec- 
essary However, unless we introduce an extra sym- 
metry the triplet r] ~ (1, 3, 0) ]20| , p^ couples also to the 
leptons. In this case we have FCNI at the tree level me- 
diated by scalars. The mass matrix is diagonalized by 
the biunitary transformation using the matrices El u: 
AP = e\{Gvs + C'v^)Er (M' = diag(TOe, m;,,mT-) and 
Vs is the vacuum expectation value (VEV) of one of neu- 
tral Higgs ( denoted cr^ in the Eq. (23) of Ref ||||) con- 



I 



and 



is the VEV of one of 



tained in the sextet |22^ 
triplets 22|) but this transformation does not diago- 
nalize Gvg nor G'vn separately. Hence, the interactions in 
Eq. (|^) are not suppressed by the lepton mass. Since Vs is 
the VEV needed only in order to give mass to the charged 
leptons it is not necessarily of the order of the hundred 
GeV but may have a value of a few GeV. Hence at this 
stage the only constraints in the Yukawa parameters G 
and G" come from perturbation theory: |G'p/47r < 1, 
|G"P/47r < 1. Since G and G' can be arbitrary complex 
matrices, the matrices defined in Eq. (||) are non-unitary 
matrices. 

We observe that in order to obtain e;^^ and appro- 
priate to solve solar or atmospheric neutrinos anomalies 
it is necessary that x,y,x',y' be not too small (that is 
not the case of grand- unified theories as discussed above) . 
The charged scalars can not be so heavy. In fact, even 
the constraints coming from the neutral kaon system do 
not impose severe restriction to the mass of these scalar: 
■mfi+ ^ nii^a ~ 150 GeV (/i+ and hP denote any of the 
singly charged and neutral scalars) if at the same time 
we take into account the mixing angles in the quark sec- 
tor 2gl. It means x ^ y ^ 0.8 already give a value 
of e^^ ~ 0.1 which is needed for solving the atmospheric 
neutrino anomaly pl] | . A fine tuning of the mixing angles 
keep e'^p almost zero. However, it is still possible that 
X, y, x' , y' be larger than one because or h'^ contribut- 
ing to the neutral kaons parameters are not necessarily 
the same component than the scalars in Eqs. (|^). The 
same analysis is valid for the anti-neutrino parameters in 
Eqs. (I). 

All scalars in Eq. (g) are still symmetry eigenstates. 
Hence, they are linear combinations of the mass eigen- 
states. It means that some of their components couple to 
quarks and induce contribution to the neutral kaon pa- 
rameters. However, those contributions depend on mix- 
ing angles in the u- and d- quarks mixing matrices so they 



do not impose constraints on the couplings in Eq. (|^). 
They will be constrained by some phenomenological pro- 
cesses as discussed later in Sec. ^ 



IV. MULTI-HIGGS DOUBLETS EXTENSION OF 
THE STANDARD MODEL 

Let us consider the standard model plus an arbi- 
trary number of scalar doublets $j = {(j)f , (j)^)'^ , i — 
1,2, ...TV > 3 In the lepton sector the Yukawa interac- 
tions are (with massless neutrinos) 

-^Y = Y1 {""L ^R^t + eIt\Er Ir 4>i) + H.c 

i 

(12) 

as in the previous model we have redefined the fields as 
in Eq. (H). The mass matrix for the charged leptons 



M' = '^i{vi/\/2)T\ is diagonahzed as follows: 



e\M^Er = m\ 



(13) 



with M = diag(r7ie, JTi/j, 'Tir). Hence, the unitary ma- 
trices El^r diagonalize M' but not each of the V\ sepa- 
rately. Although we have redefined the neutrino fields in 
the charged currents coupled to the vector bosons W , the 
same is not possible in the interactions with . Hence, 
even with massless neutrinos we cannot avoid, in gen- 
eral, to have mixing in the charged current mediated by 
scalars in the lepton sector as well. Moreover, even if 
neutrinos are massless, there are flavor changing interac- 
tions mediated by neutral scalars in the charged lepton 
sector. 

The charged currents in Eq. ( |l^ ) can be written in 
terms of the physical charged scalar TLf (defined as (j)^ — 



H.c. 



where 



kOL 



(14) 



(15) 



where a = e, /z, t and k = 1,2,3. The matrices Vj are not 
unitary matrices since F' are in general arbitrary complex 
matrices, with |rp/47r < 1. 
In this case we have 



Explicitly 



E 



(vj)^^(v.).,-(vt)^^(v,)^, 



(18) 



with = A/^ /M^ . . Note again that in order to obtain 
the relevant FCNI parameters for solar neutrino analysis 
it is sufficient to change appropriately the indexes /i — > 
e in the Eqs. ( p7| ) and (p^). In this case eap = —^ap 
(similarly for the e's) and the constraints coming from 
/i — > 67 and other processes are similar as for the previous 
model as it will be shown in the next section. 



V. CONSTRAINTS ON FCNI PARAMETERS 

In order to be sure that a realization of the oscilla- 
tions induced by FCNI can happen in a quantitative way 
and a compatibility of data and related predictions be 
achieved we have finally to investigate the constraints 
to these FCNI that arise from the charged lepton decay 
measurements and from the absence of violation of the 
lepton number conservation laws | p7[ |. 

The /i-decay is measured with large precision and it is 
consistent with the standard electroweak model predic- 
tions. This imposes severe constraints on exotic inter- 
actions involving neutrinos. Considering the 331 model, 
H2 and mediate processes which should be summed 
coherently to the standard model contribution. The in- 
terference contribution of iJ^ and 77^ with the standard 
model (SM) prediction, for the the left-handed /z-decay 
measurement, leads to |28||: 



Re 



:2(^Ll)ip V' + (A:[i)e2(/CLL)i^ x' 

< (Gf /y2)M2,A^, 



(19) 



where is the combined experimental and theoretical 
error of /^-decay. Following the recipe of Ref. ||2^ we 
compute the Gf.*^ value 



r^SM 



40) 



V2M^ (1 - M^/Aq) (1 - Ar) 



(20) 



In Eq. (^0|) the fine structure constant a(0) is very ac- 
curately known to be 1/137.036 and Ar is the radiative 
SM correction and is found to be Ar = .0349 T -0019 ± 
.0007 |2^]. However, in computing the largest error 
comes from Mw measurement with current value given 
as Mw — 80.39 ± .06 GeV |Q and a much smaller error 
comes from Mz = 91.1867 ± .0020 GeV 

Using this and propagating the error of Gf.^^ in the 
error of /z-decay width we get that at la, A^ = 8.6 x lO^'^. 
The mainly error came from the Gf.*^ uncertainty. Also 
in the case of r-decay we need to take into account the 
T-lifetime measurement error and r-mass measurement 
error, then A^ = 9.5 x 10~^ at la. 
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We can consider also the r-decay measurement 



Re 



(^LL)e3(^LL)lr v' + (/C[^)e3 (/Cll)i^ x' 



(21) 



which constraint is obtained from Eq. (|T^) substituting 
appropriately the sub-indexes and imposing A,- ~ 9.5 x 
10"'^ which leads to less stringent limits. 

We can rewrite, using the symmetric (anti-symmetric 
) properties of ICll {^'ll )> 



Re 



and 



Re 



-{^'IMic'^lW y' + (4l)2.(^ll)im ^' 

-(^LL)3e(^LL)lr V' + (^lL)3e(/CLL)lr x' 



(22) 



< {Gf/V2)MI,/^, 



(23) 



These constraints can give limits on the contribution of 
H2 and 77^ to the e^^ defined in Eqs. (|5b|) . Respectively 
the ^ decay and the r decay can give limits on and 



£33. Nevertheless, Eqs. ( |22| ) and (^3|) can easily exhibits 
cancellations among the contributions which can elimi- 
nate the constraints. 

There are also other contributions of the charged- 
lepton decay mediated by H2 and which do not sum 
coherently to the standard electroweak process 



(24) 



where j = /i, r and i ^ a and j ^ (3. These relevant 
constraint is to the e^^ defined in Eq. (5b) is given for 



1 = 1 and a = 1 in Eq. (^4|). Then the contribution 

to e^p is 



-1/3 1 ff,- 



(25) 



where /3 = 2, 3. The contribution of 77^ vanishes when 
a = 1 because the anti-symmetric character of ICll- 

Also we have the contributions of H2 and that are 
constrained 



(26) 



In the 331 model there are contributions to the right- 
handed charged-lepton decays mediated by V^, and 
?7f do not sum coherently to the standard process. The 
relevant constraint is 



Eo,/3 |(/Ct)„,(/C),0|^'^ < (Gy2)M^A^ 



(27) 



with a and /3 are related to the final neutrino states, j and 
i is related to the initial and final leptons respectively. If 
i = 3 [j = 2) implies r (^)-decay. Recall that although /C 
is a unitary matrix, no symmetry relation exists in K-lr 
or /C^^ couplings. Each of elements of the sum defined 
in Eq. ( p7| ) also need to be constrained to be lower then 
the right side. 

For example, the combination of product of matrices 
K, given in Eq. (|27| did not appear in the definition of 
(Vide Eq. (^b|)), but some of matrices element can 
be constrained to be small. The element {JO)aj with 
j > 1, or (/C)jQ did not appear in the Eq. (^), but can 
be constrained in the r — > /i7 decay (see below). Similar 
arguments apply for the K-lb. or /C^^ couplings. 

Concerning the contributions to the charged lepton 
decays arising in the multi-Higgs models of Section IV, 
there are no contributions to be coherently summed to 
the standard model contribution to the charged lepton 
decays. Therefore the relevant constraints coming from 
the charged lepton decay measurements can be written 
as 



< {M^Gf/V2YA 



(28) 



with a and [3 related to the final neutrino states, k and 
I is related to the initial and final leptons respectively. 
Also each term of sum of Eq. (28) is constrained to be 
smaller then the right side of Eq. ( p8|) . 

Now, from Eq. ([l6|) we can rewrite 



1 



1/3 



Then comparing Eq. (g9|) and (|2^) we get 

I ^ |2 



< Aa/64 



(29) 



(30) 



where a = 2, 3 and any /3. 

These is a stronger constraint on the multi-Higgs 
model. Only the element e\i is not constrained. 

Lepton number conservation tests pTf can generate 
strong limits to FCNI parameters. In 331 model and 
in the multi-Higgs model there are constraints coming 
from ^ e7, r — > 67, r — > /i7 and ^ — > eee. In 
the multi-Higgs doublets model there are contribution 
to the /X ^ 3e decay mediated by the neutral Higgs fla- 
vor changing interactions given in Eq. (^2|) it constrains 
the same matrix elements that appear in the charged cur- 
rents in Eq. (^2|) so, the mass of the neutral Higgs bosons 
and the respective mixing angles. 

In 331, appear contributions to the /i — *■ 67 decay via 
the exchange of a gauge vector boson V {^ji ^ Vv^ 
67); via exchange of charged Higgs iJf and rj^ {fi 
Hi{rii)va 67), as well of iJ^ and 77^ ; via exchange 
of doubly charged gauge boson U~~ (the interactions are 
not showed here) |^ and via exchange of doubly charged 
Higgs Ht+ {H++). 



The amplitude of the contribution of gauge vector bo- 
son V is proportional to ^2aA^li and vanish because 
/C is a unitary matrix. 

The contribution of charged Higgs and r/f can be 
written as 



< (487r/a(0))M4,Bi?(^ ^ 67) 



(31) 



where BR{^ ej) < 1.2 x 10"^^ Q is the branch- 
ing ratio of the lepton violating process — s- ej. Note, 
however, that the coupling {K,LR)a2 and (/C^fl)c(2 do not 
appear in Eqs. (|5^) and (^) and therefore the limit above 
does not imply direct constraint on e^^ and e^^. 

The same happens with the contributions of H2 and 
■ Only the sum of some matrix elements is constrained: 



< {A8TT/aiO))M^BR{^i -> 67). 



(32) 



Even considering that the matrix ICll {^'ll) sym- 
metric (antisymmetric), their combination appearing in 
Eq. (^) does not contribute to the definition of e|^^. Fur- 
thermore the constraint involves a coherent sum in which 
we can find cancellations making it weaker. 

Finally there is a contribution of the doubly charged 
vector boson U , /i U — s- 67; since the mass 
of the charged leptons are quite different, the amplitude 
only constrains | Eq ^A'a^ae'^Q/^'^lf-- vf\th nia = 
me,m^,mT- Since the m,- dominates, the decay con- 
strains only ICjj^^lCse- These contributions should be co- 
herent summed to the and to the . The former 
have the couplings ICrr (See Eq. (||)) that did not appear 
in the Eqs (5a) and (|5b|) and the later Kll that appear in 
those equations, but instead of H2 mass in the denomi- 



nator of the Eqs (5a) and (5b) you have the H2 mass 



in the denominator. You can choose to put the masses 
of H2" field large enough to suppress the constraints on 

The /i ^ 3e decay arise by the interaction with the 
doubly charged Higgs H^^ (^^2"^) ^'^^ t)y the interac- 
tion with the doubly charged vector boson C/++. The 
Higgses H^^ {H^^) constrain the matrix JCll {K-rr) 
and doubly charged vector boson U'^^ which constrains 
the K, mixing matrix. The matrix elements involved are 



K/^i and /Cei and only the first one appears in Eq. (|7aj) 



Note, however, that the doubly charged vector and scalar 
bosons do not contribute to e^^. Any constraint can 
be evaded by imposing that their masses are sufficiently 
large. 



At first sight the process yLv,, 



would impose 



constraints on {K,lr)i^ and {Ki'i^j^^^ and other param- 
eters involving the ^. However, this process has a cross 
section which is 0.90 ± 0.20 times the prediction of the 
V - A theory || (or 0.98 ± 0.18 H) which does not 



imply therefore strong constraints on the related cou- 
plings. The doubly charged scalars 11^2 may contribute 
to the muonium to anti-muonium conversion with the 
same strength as the doubly charged vector bilepton. It 
depends on the values of the parameters such as vacuum 
expectation values of the model |Q . 

Finally we mention that the FCNI parameters can be 
constrained using only charged leptons decays and some 
violating lepton number processes. We get that the con- 
tributions of charged Higgs scalars H2 and 77^ to ef^, 
with /3 = 2, 3 and £^3 are severely limited. Also some 
combinations of couplings are constrained to be very 
small like Eqs (||), dH), (|l]) and (js^) although have 
no direct comparison with the relevant couplings of e^^ 
and e^^. For the other couplings like K-lr, T^'lr and K. 
we have a compromise between different processes but no 
direct limit can be obtained. 

For the case of multi-Higgs, the situation is different. 
Only the charged lepton decays impose a stronger limits 
on the e^^ and e^^. Unless you enlarge the model, e. 
g. adding a set of singly charged Higgs which can made 
less drastic the constraints, we have no hope to get FCNI 
parameters large enough to get a relevant oscillation. 

Also is worth to mention that the e^^ and e^^ pa- 
rameters are essentially different, then even in the case 
of massless neutrinos, the neutrinos and anti-neutrinos 
have different evolution in the matter. This is not the 
case in the i?-parity violating supersymmetric models. 



VI. CONCLUSIONS 

Both models discussed above can generate FCNI with 
strictly massless neutrinos with strength sufficient to in- 
duce appreciable effects both for solar as well as atmo- 
spheric neutrinos. These effects can either be confirmed 
or ruled out as the collected data sample increases and 
the solar and atmospheric results become more precise. 

If confirmed, it is necessary to know what is the actual 
theoretical framework that can account for these type of 
FCNI. In fact, there are very few models in the literature 
that introduced FCNI while keeping neutrinos massless. 
If, on the other hand, future data points undoubtedly in 
the direction of conventional neutrino oscillation (in vac- 
uum or resonantly enhanced by matter), FCNI present in 
these models will have to be suppressed and these models 
will either be ruled out or their FCNI couplings severely 
restricted by data. 

It is worth to mention that a very interesting interme- 



_ f 

diate scenario occurs when oscillations induced by e 



and co-exist with conventional neutrino oscillations. 
It was pointed out that in this situation simultaneously 
solution to the solar, atmospheric and LSND observa- 
tions l^sl can be obtained in a context of only three neu- 
trino families p6|], dispensing therefore the introduction 



Q/3 



7 



of a fourth light (electroweak singlet) neutrino. 

Constraints from several processes previously analyzed 
can generate limits to the FCNI parameters e^^ and e^^. 
We see above that the 331 model survives like a model 
able to generate FCNI parameters compatible with the 
constraints and a viable model of FCNI and still keeping 
the neutrinos massless. This is a different situation from 
what happens in i?-parity broken supersymmetric models 
where the combination of matrices entering in the defini- 
tion of the FCNI parameters appear also contributing to 
constrained lepton violating processes. It is also different 
from what happens in the Zee model. There the unique 
anti-symmetric coupling makes it impossible to avoid the 
constraints. 

So far, neutrinos remain massless. In both models we 
can get massive neutrinos simply by adding the right- 
handed components and the neutrino masses are arbi- 
trary. In the 331 model we have still the possibility that 
one of the neutral component of the scalar sextet |^ ( 
denoted in the Eq. (23) of Ref |Q) get a non- vanish 
VEV, here denoted by vi, giving to the neutrinos a ma- 
jor ana mass term 



Gvi 



(33) 



where is a possible Majorana mass term for the right- 
handed singlets. The neutrino masses are still arbitrary 
but there are FCNI in the scalar sector. It means that 
there are Yukawa couplings that are not proportional to 
the neutrino masses. In both models we can also add 
several neutral or doubly charged scalar singlets only to 
get calculable neutrino masses [ p7| , ^ . Or, it is possible 
to get calculable neutrino masses breaking explicitly the 
total lepton number in the scalar potential In all 

these possibilities the neutrino masses do not depend on 
the parameters entering in the e and e'. One example 
of new contribution to the FCNI effect, when the right- 
handed neutrino is introduced, is showed in Fig. H. 



V 



'(3R 



FIG. 4. Contribution if right-handed neutrinos are added. 

Hence, we see that there are no strong constraints on 
all the parameters of the 331 model. And from this we 
hope that show a viable model for flavor changing neu- 
trino interactions with strictly massless neutrinos that 
survive the constraints from charged lepton decays and 
some violating processes like /i — > 67, ^ eee. 
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